AFIT/GE/ENG/99M-2 1 


LOW  GRAZING  ANGLE 
TERRAIN  BACKSCATTERING 
COEFFICIENT  (a0)  GENERATOR 

THESIS 


William  N.  O’Connor 
Captain,  USAF 
AFIT/GE/ENG/99M-2 1 


Approved  for  public  release;  distribution  unlimited 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  121 5  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

March  1999  Master’s  Thesis 

4.  TITLE  AND  SUBTITLE 

LOW  GRAZING  ANGLE  TERRAIN  BACKSCATTERING  COEFFICIENT 
GENERATOR 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

William  N.  O'Connor,  Captain,  USAF 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Institute  of  Technology, 

2950  P  Street 

WPAFB,  OH  45433-7765 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/GE/ENG/99M-21 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFRL/SNAS 

Attn:  Capt  Kelce  Wilson 

2010  Fifth  Street,  Bldg  23,  Area  B 

Wright-Patterson  AFB,  OH  45433 

19371  255-6329/DSN  785-6329 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

Dr.  Vittal  Pyati  /  255-3636  x4576,  DSN  785-3636  x4576  /  vittal.pyati@afit.af.mil 

12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

12b.  DISTRIBUTION  CODE 

1 3.  ABSTRACT  (Maximum  200  words) 

A  near-horizon  probabilistic  terrain  backscatter  coefficient  generator  is  introduced,  combining  modified  sea-state  curves  with 
Weibull-type  probability  density  functions.  A  total  of  nine  terrain  type  classifications  (soils  &  rocks,  trees,  grasses,  shrubs, 
short  vegetation,  roads,  urban  areas,  dry  snow,  and  wet  snow)  are  supported,  for  seven  radar  bands  (L,  S,  C,  X,  Ku,  Ka, 
and  W),  and  both  linear  (HH,  W)  and  cross  (HV,  VH)  polarizations.  The  Weibull  functions  are  based  on  measured  sigma 
and  match  the  published  data.  The  functions  also  provide  a  probabilistic  nature  to  the  generator.  Modified  sea  state  curves 
are  used  in  place  of  the  Weibull  functions  for  near-horizon  incidence  angles  due  to  the  dearth  of  measured  data  at  this  range. 
The  sea  state  approach  is  based  on  the  theory  that  for  sigma,  similar  surface  roughness  of  two  different  terrain  types  results  in 
similar  sigma  curve  shapes.  Generator  testing  and  validation  used  available  measured  sigma  for  three  terrain  types.  The 
testing  did  not  reject  the  modified  sea  curve  usage  and  allowed  system  implementation  for  the  remaining  terrain  types. 

14.  SUBJECT  TERMS 

Radar,  Radar  Backscatter,  Low  Grazing  Angle,  Ground  Clutter,  Clutter  Suppression 

15.  NUMBER  OF  PAGES 

86 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION 
OF  REPORT  OF  THIS  PAGE 

Unclassified  Unclassified 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

UL 

Standard  Form  298  {Rev.  2-89)  (EG) 

Prescribed  bv  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 


AFIT/GE/ENG/99M-2 1 


LOW  GRAZING  ANGLE  TERRAIN 
BACKSCATTERING  COEFFICIENT  (a°)  GENERATOR 

William  N.  O’Connor,  B.S. 

Captain,  USAF 


Approved: 


(P-  _  3  Mg*  7  f 

Dr.  Vittal  P.  Pyati  Date 

Thesis  Advisor 


AFIT/GE/ENG/99M-2 1 


LOW  GRAZING  ANGLE  TERRAIN 
BACKSCATTERING  COEFFICIENT  (a0)  GENERATOR 

THESIS 

Presented  to  the  Faculty  of  the  Graduate  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Command 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  in  Electrical  Engineering 

William  N.  O’Connor,  B.S. 

Captain,  USAF 

March  1999 


Approved  for  public  release,  distribution  unlimited 


Acknowledgments 


I  would  like  to  thank  my  thesis  advisor,  Dr.  Vittal  Pyati,  and  my  thesis  committee, 
Maj  Michael  Temple  and  Dr.  Andrew  Terzuoli,  for  their  guidance  and  assistance  during 
this  project.  They  have  shown  great  patience  and  foresight,  and  many  times  they  have  kept 
me  from  going  down  a  wrong  path. 

My  sponsor,  Capt  Kelce  Wilson,  was  no  less  important  to  the  success  of  this 
research.  His  knowledge  and  expertise  on  this  subject  helped  greatly  and  I  appreciate  all 
the  time  he  took  from  his  busy  schedule  to  answer  questions  or  to  listen  to  crazy  ideas. 

The  guys  of  the  Communications  Lab,  Rick,  Andy,  Andy,  Perry,  Todd,  John,  Dan, 
Patrick,  Dave,  and  Mike.  Together,  we  spent  many  a  sleepless  night  working  on  papers 
and  homework  while  always  keeping  alive  a  sense  of  humor.  Even  with  all  the  hard  work 
and  pain  we  went  through,  I  will  miss  these  times. 

Finally,  to  my  family,  Neil,  Patricia,  Patrick,  Michael,  and  the  pets;  Thank  you  for 
being  the  support  I  could  always  depend  on.  Words  can  not  fully  describe  what  you  have 
done  for  me,  so  I  will  simply  say,  Thank  you. 


u 


Table  of  Contents 


Acknowledgments . ii 

List  of  Figures . . 

List  of  Tables . . . vii 

Preface . viii 

Abstract . ix 

l.  Introduction . 1 

1.1  Background . 1 

1 .2  Problem  Statement . 2 

1 .3  Summary  of  Current  Knowledge . 2 

1 .4  Thesis  Overview . 4 

n.  Theory . 6 

2.1  Overview . 6 

2.2  Calculation  of  Backscatter  Coefficient  (a0) . 8 

2.3  Terrain  Dependence  of  g° . 1 1 

2.4  Simulated  Surface  Model . 12 

2.5  Analysis  of  Available  Grazing  Data . 29 

2.6  Currently  Available  Low-Grazing  Terrain  Backscatter  Models . 30 

2.7  Summary . 33 

m.  Methodology . 34 

3.1  Introduction . 34 

3.2  Sea  Model  Theory . 35 

3.3  Wilson/Mediavilla  Model . 40 

3.4  Model  Integration . 40 

3.5  Model  Validation . 44 

3.6  Summary . 47 

IV.  Results  and  Discussion . 48 

4.1  Overview . 48 

4.2  Sea  State  decisions . 48 

4.3  Model  Results  for  selected  tests . 49 

V.  Conclusions  and  Recommendations . 64 


iii 


Appendix  A.  Terrain  Backscatter  Code  in  MATLAB® . 66 

A.  1  Introduction . 66 

A.2  WMOC  Function . 66 

A.3  SEAMODEL  Function . 69 

Bibliography . 72 


List  of  Figures 


Figure  1. 
Figure  2. 
Figure  3. 
Figure  4. 


Figure  20. 
Figure  21. 


Grazing  angle  vs.  Incidence  Angle . 2 

Radar  at  Horizon  receiving  a  distant  object . 4 

Illumination  of  an  area . 10 

Random  Surface . 13 

Figure  5.  Scattering  off  of  Random  Surface  [14] . 14 

Figure  6.  Visible  Facets  in  dB  (Grass) . 15 

Figure  7.  Average  Visible  Facet  Slope  Angle  (Grass) . 16 

Figure  8.  Vertical  Coefficient  (Grass) . 16 

Figure  9.  Horizontal  Coefficient  (Grass) . 17 

Figure  10.  Total  Coefficient  in  dB  (Grass) . 18 

Figure  11.  Visible  Facets  in  dB  (Trees) . 20 

Figure  12.  Average  Visible  Facet  Slope  Angle  (Trees) . 20 

Figure  13.  Vertical  Coefficient  (Trees) . 21 

Figure  14.  Horizontal  Coefficient  (Trees) . 21 

Figure  15.  Total  Coefficient  in  dB  (Trees) . 22 

Figure  16.  Sea  Waves . 23 

Figure  17.  Visible  Facets  in  dB  (Sea) . 24 

Figure  18.  Average  Visible  Facet  Slope  Angle  (Sea) . 24 

Figure  19.  Vertical  Coefficient  (Sea) . 25 

Horizontal  Coefficient  (Sea) . 25 

Total  Coefficient  in  dB  (Sea) . 26 

Figure  22.  Measured  values  of  W  and  HV  coefficients  for  ks  <  2  [7] . 28 

Figure  23.  Measured  values  of  W  and  HV  coefficients  for  ks  >  5  [7] . 28 

Figure  24.  Georgia  Tech  Land  Model  vs.  Sea  Clutter  Model, . 32 

Figure  25  Sea  Model  Curves  for  5  Different  Sea  States . 39 

Figure  26.  Trees,  Ku-Band,  Before  averaging  to  plateau . 42 

Figure  27.  Trees,  Ku-Band,  After  averaging  to  plateau . 42 

Figure  28.  Grasses,  S-Band,  Before  averaging  to  plateau . 43 

Figure  29.  Grasses,  S-Band,  After  averaging  to  plateau . 43 

Figure  30.  Comparison  of  the  final  model  with  ARL  data  (Trees  Coniferous,  W-Band,  W- 

polarization) . 45 

Figure  31.  Comparison  of  the  final  model  with  ARL  data  (Trees  Deciduous,  W-band,  W- 

polarization) . 45 

Figure  32.  Comparison  of  the  final  model  with  ARL  data  (Grasses,  W-Band,  VV-polarization) ....  46 
Figure  33.  Comparison  of  the  final  model  with  ARL  data  (Short  Vegetation,  W-band,  VV- 

polarization) . 46 

Figure  34.  L-Band . 50 

Figure  35.  Ka-band . 51 

Figure  36.  Ku-Band . 52 

Figure  37.  X-Band . 52 

Figure  38.  C-Band . 53 

Figure  39.  Ku-Band . 54 

Figure  40.  C-Band . 55 

Figure  41.  S-Band . 55 

Figure  42.  Ku-Band . 56 


Figure  42.  Ku-Band . 56 

Figure  43.  L-Band . 57 

Figure  44.  Ka-Band . 58 

Figure  45.  X-Band . 58 

Figure  46.  L-Band . 59 

Figure  47.  X-Band . 60 

Figure  48.  C-Band . 61 

Figure  49.  W-Band . 61 

Figure  50.  S-Band . 62 

Figure  51.  Ka-Band . 63 


\ 


vi 


List  of  Tables 


Table  1.  A%  Change  in  factors  (Grass) . 19 

Table  2.  A%  Change  in  factors  (Trees) . 22 

Table  3.  A  %  Change  in  factors  (Sea) . 26 

Table  4.  ARL  measured  data,  Mean  c°  and  Standard  Deviation  in  dB . 30 

Table  5.  Input  parameters  for  Sea  Clutter  Model . 35 

Table  6.  Internal  variables  for  Sea  Clutter  Model  and  dependence . 35 

Table  7.  Sea-State  Surface  Descriptors . 36 


vii 


Preface 


This  research  develops  a  near-horizon/low-grazing  angle  backscattering  coefficient 
(a0)  generator  for  several  radar  bands,  types  of  terrain,  and  wave  polarizations. 

The  low-grazing  a0  generator  is  a  fusion  of  two  separate  models:  a  low-incidence 
angle  g°  generator  that  returns  the  probabilistic  means  of  g°  for  a  given  terrain,  frequency, 
and  incidence  angle;  and  a  low-grazing  sea  surface  backscatter  model.  The  theory  that 
low-grazing  angle  G°  depends  more  on  surface  roughness  than  on  surface  composition 
allows  for  this  sea  model  usage. 

The  generating  function  is  implemented  in  MATLAB®,  and  is  updateable  when 
new  low-grazing  data  becomes  available.  The  function  is  available  in  FORTRAN  from  the 
author. 
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Abstract 

A  near-horizon  probabilistic  terrain  backscatter  coefficient  (a0)  generator  is 
introduced,  combining  modified  sea-state  curves  with  Weibull-type  probability  density 
functions.  A  total  of  nine  terrain  type  classifications  (soils  &  rocks,  trees,  grasses,  shrubs, 
short  vegetation,  roads,  urban  areas,  dry  snow,  and  wet  snow)  are  supported,  for  seven 
radar  bands  (L,  S,  C,  X,  Ku,  Ka,  and  W),  and  both  linear  (HH,  VV)  and  cross  (HV,  VH) 
polarizations.  The  Weibull  functions  are  based  on  measured  G°  and  match  the  published 
data.  The  functions  also  provide  a  probabilistic  nature  to  the  generator.  Modified  sea  state 
curves  are  used  in  place  of  the  Weibull  functions  for  near-horizon  incidence  angles 
(75°  -  90°)  due  to  the  dearth  of  measured  data  at  this  range.  The  sea  state  approach  is  based 
on  the  theory  that  for  G°,  similar  surface  roughness  of  two  different  terrain  types  results  in 
similar  G°  curve  shapes.  Generator  testing  and  validation  used  available  measured  G°  at  85° 
for  three  terrain  types.  The  testing  did  not  reject  the  modified  sea  curve  usage  and  allowed 
system  implementation  for  the  remaining  terrain  types. 
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LOW  GRAZING  ANGLE  TERRAIN 


BACKSCATTERING  COEFFICIENT  (a0)  GENERATOR 

I.  Introduction 

1.1  Background 

The  backscattering  coefficient,  a0,  is  a  measure  of  the  electromagnetic  energy 
scattered  by  a  distributed  target  in  the  receiver  direction.  In  general,  o°  depends  upon  the 
target's  shape,  dielectric  constant,  viewing  geometry,  and  the  attributes  of  the  incident  and 
scattered  electromagnetic  waves  (wavelength,  polarization,  and  incidence  angle)  [1]. 

For  terrain,  the  value  of  g°  depends  upon  its  features,  such  as  water  content, 
composition,  degree  of  roughness,  and  the  location  of  scattering  centers.  Generally,  terrain 
scattering  centers  are  randomly  distributed  and  vary  in  position  with  frequency  and 
incidence  angle,  giving  a0  its  random  quality  [1].  The  variation  of  the  backscattering 
coefficient  justifies  its  description  as  a  mean  value  calculated  from  several  measurements  at 
a  specific  incidence  angle,  frequency,  and  wave  polarization.  The  associated  standard 
deviation  of  the  measurements  may  be  calculated  and  used  as  well. 

For  low-grazing  angles  near  the  horizon,  the  data  available  is  very  scant.  Because 
of  this  lack  of  data,  deterministic  models  have  been  developed  based  on  small  subsets  of 
the  available  data  that  would  fit  the  determined  model. 

The  terms  “low-grazing  angle”  and  “high-incidence  angle”  are  used  interchange¬ 
ably.  This  exchange  in  terms  is  owed  to  the  popularity  of  the  “low-grazing”  term  in  the 
scientific  community  and  in  early  G°  research  while  the  term  “high-incidence”  is  used  by 
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the  two  primary  sources  of  data.  In  this  work,  for  continuity  with  past  studies,  “low- 
grazing”  will  describe  the  region  of  research  and  the  measured  data.  However,  all  angles 
described  use  incidence  nomenclature  for  continuity  with  the  two  data  sources  (University 
of  Michigan,  RADLAB  and  Army  Research  Laboratory,  ARL) 


0  HORIZON 
20 


0  90 

NADIR  NADIR 


Figure  1.  Grazing  angle  vs.  Incidence  Angle 


1.2  Problem  Statement 

Develop  a  low-grazing/near-horizon  c°  generator  for  several  terrain  types, 
frequency  bands,  and  wave  polarizations.  This  generator  will  operate  in  conjunction  with 
the  backscatter  coefficient  generator  [2]  developed  for  incidence  angles  0°  <  0  <  80°.  The 
generator  will  output  data  for  the  nine  terrain  types  specified  by  RADLAB  [1],  seven  radar 
bands,  and  three  wave  polarizations.  When  new  low-grazing  angle  data  becomes  available, 
the  generator  must  be  updateable  and  generate  G0,s  with  the  same  probabilistic  distributions 
as  the  measured  data. 


1.3  Summary  of  Current  Knowledge 

Very  little  measured  data  is  available  for  low-grazing  angle  G°  measurements.  Vast 
majority  of  studies  have  created  deterministic  models  based  only  on  the  limited  amount  of 
data  they  measured  [3,  4,  5],  Army  Research  Laboratory  (ARL)  measured  W-band 
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backscatter  coefficients  from  85°  <  0  <  89°  for  European-type  terrain  (trees,  grass  fields, 
dirt,  and  snow)  [6].  Unfortunately,  the  measurement  sample  sizes  and  number  of  terrain 
types  measured  were  limited.  RADLAB  has  been  compiling  a0,s  from  numerous  studies 
and  creating  a  large  o°  database  complete  with  means  and  standard  deviations  for 
numerous  terrain  types  and  incidence  angles.  The  RADLAB  data  compilation  is  the  basis 
for  the  lower  incidence-angle  c°  generator  being  concurrently  developed  [2].  However, 
there  is  no  data  available  from  RADLAB  for  incidence  angles  0  >  80°. 

1.3.1  Beckmann  Method. 

In  general,  when  low-grazing  data  was  limited  or  unavailable,  almost  all  prior 
models  use  a  squared  sinusoid  to  fill  in  the  data  near  the  horizon. 

<7°(0)  =  (To(o)-sin20  (for  grazing  angles)  (1) 

<7o(0)  =  (T0(0)cos26  (for  incidence  angles)  (2) 

where  G°(0)  is  the  terrain  a0  value  at  nadir. 

While  this  method  has  a  good  fit  of  measured  data  for  a  limited  number  of  terrain 
types  and  incidence  angles,  it  has  been  rejected  for  its  behavior  at  the  horizon  [3,  4,  6], 
When  0  is  at  the  horizon  (0  =  0°  for  grazing,  90°  for  incidence),  a0  equals  zero.  A  number 
of  studies  that  have  measured  a0  at  the  horizon  have  detected  a  small  amount  of  return 
energy  [3, 4, 7],  This  return  energy  usually  is  due  to  a  distant  object,  or  objects,  in  view  of 
the  receiver  (a  tree,  long  grass,  sea  waves,  or  a  distant  building).  Figure  2  demonstrates  this 
condition. 
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Figure  2.  Radar  at  Horizon  receiving  a  distant 
object 


1.4  Thesis  Overview 

A  backscatter  generator  for  incidence  angles  0°  <  0  <  80°,  based  on  the  probabilistic 
distribution  of  measured  a°'s  ,  is  available  [2]  (to  be  referred  as  the  Wilson/Mediavilla 
Generator  or  Gwm  GEN).  By  theorizing  that  at  low-grazing  angles,  c°  is  more  affected  by 
surface  roughness,  a  widely  accepted  sea  backscatter  model  is  used  to  generate  backscatter 
curves  for  incidence  angles  up  to,  and  including,  horizon.  The  overall  value  or  “height”  of 
the  curve  is  set  by  Gwm  GEN  and  the  sea  model  will  extrapolate  the  rest  of  the  desired  G°’s. 
[2] 

Terrain  roughness  is  specified  by  sea  wave  height.  The  roughness-wave  height 
correlation  is  done  using  available  measured  data  as  a  validating  target.  In  situations  where 
data  is  unavailable,  the  wave  height  will  be  set  by  either  the  wave  height  given  for  the 
specified  terrain  at  a  higher  frequency  (usually  W-band  since  the  only  available  data  is  in 
W-band)  or  by  using  a  Rayleigh  roughness  model.  [5] 

Chapter  2  reviews  the  theory  behind  the  measurement  of  o°  and  a  method  to 
theoretically  measure  the  effect  on  low-grazing  angle  a0  of  surface  roughness  vs  surface 


composition.  The  chapter  concludes  with  an  analysis  of  the  ARL  data  and  shows  how  it 
relates  to  the  RADLAB  data  used  in  (Twm  GEN. 

Chapter  3  will  introduce  the  sea  backscatter  model  and  validation  of  the  model 
using  Gwm  GEN  and  the  ARL  data. 

Selected  G°  generations  for  the  near-horizon  incidence  angle  range  (60°  <  0  <  90°) 
will  be  graphed  in  chapter  4. 

Conclusions  and  recommendations  for  future  study  will  sum  up  chapter  5  and  this 
research. 
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II.  Theory 


2.1  Overview 

Since  the  beginnings  of  radar,  the  detection  of  unwanted  objects,  or  “clutter”,  has 
interfered  with  the  measurement  of  desired  targets.  The  problem  is  that  there  is  no  one  type 
of  clutter.  An  airport  search  radar  tracking  aircraft  would  call  weather  returns  clutter,  while 
a  weather  radar  would  want  to  remove  the  clutter  caused  by  aircraft.  Furthermore,  a  space 
surveillance  system  would  consider  aircraft  and  weather  returns  to  be  unwanted  clutter. 
Thus,  this  gives  rise  to  the  saying:  One  radar’s  clutter  is  another  radar’s  target  [8]. 

In  this  research,  our  target  is  ground  clutter,  or  now  known  as  terrain  backscatter 
(a0).  Backscatter  is  transmitted  energy  reflected  by  ground,  or  surface,  elements  back  to 
the  radar  receiver  [9].  The  amount  of  energy  received  varies  with  the  terrain  type  and  the 
angle  the  transmitted  beam  hits  the  surface.  On  most  radar  systems,  backscatter  is  just  an 
annoyance  to  be  worked  around.  However,  in  some  applications,  mostly  military, 
backscatter  can  mask  targets,  making  the  targets  more  difficult  to  detect  [10]. 

To  counteract  the  backscatter  effects  on  ground-based  radar  systems,  one  of  two 
methods  is  usually  applied.  The  first  method  uses  zero-Doppler  effects  to  identify  and 
remove  clutter  effects  [9].  This  means  that  the  system  removes  any  return  that  does  not 
appear  to  be  moving.  While  this  method  is  effective  in  removing  clutter,  it  has  been  found 
that  slow  moving  targets,  in  ground  applications,  are  also  removed,  degrading  the 
effectiveness  of  the  system  [10]. 

A  new  method  that  has  started  to  take  hold  is  predicting  the  backscatter  for  a  given 
terrain  surface  at  a  given  angle.  Because  this  method  requires  previous  knowledge  of 
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backscatter  energy  from  specific  terrain  types  as  well  as  knowledge  of  the  terrain  surface 
that  will  be  in  the  radar  bore-site,  this  method  has  been  slow  in  gaining  acceptance. 
However,  with  the  availability  of  DTED  (Digital  Terrain  Elevation  Data)  maps  for  the 
world  as  well  as  the  development  of  backscatter  models,  this  method  is  becoming  more  and 
more  prevalent  [5, 10, 1 1]. 

This  method  is  not  without  its  problems.  Because  of  the  number  of  variables  that 
can  affect  backscatter  (terrain  type,  surface  roughness,  frequency,  and  moisture  content  to 
name  a  few),  the  actual  energy  return  can  vary  from  measurement  to  measurement  [10,  1 1]. 
Usually,  the  amount  of  energy  recorded  is  a  mean  of  the  measurements.  Unfortunately,  the 
models  that  are  developed  from  the  recorded  data  are  usually  deterministic  (ie.  for  a  given 
angle,  terrain  type,  frequency,  etc. . .  the  same  result  is  returned).  The  lack  of  a  probabilistic 
nature  in  the  model  does  not  reflect  the  actual  measured  return  energy  [2,  10]. 

As  more  and  more  data  is  collected  and  recorded,  probabilistic  models  have  started 
to  develop.  These  models  add  the  variance  that  exists  in  backscatter  energy.  Most  notably, 
the  data  collection  from  RADLAB  of  the  University  of  Michigan  has  resulted  in  not  fewer 
than  3  probabilistic  models  (RADLAB,  Wilson,  &  Wilson/Mediavilla)  [1,  2,  10]. 
Unfortunately,  this  data  collection  only  exists  in  the  incidence  angle  range  from 
0°  <  0  <  80°.  The  Army  Research  Laboratory  (ARL)  measured  backscatter  energy  in  the 
"shadow"  or  grazing  region  (80°  <  0  <  90°)  for  a  small  number  of  terrain  types  at  95  GHz 
[6]. 

This  research  will  develop  a  probabilistic  model  for  the  grazing  region  by 
combining  an  available  grazing  angle  deterministic  model  with  one  of  the  above  mentioned 
probabilistic  models  (Wilson/Mediavilla)  [2]. 
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In  this  chapter,  a  sample  G°  will  be  calculated  from  the  radar  range  equation  and  the 
effect  of  terrain  on  o°  will  be  documented.  Then,  three  sample  surfaces,  two  random  sets 
of  grass  or  trees,  and  one  sea  waves  set,  of  varying  heights  will  demonstrate  geometric  and 
reflectivity  effects  on  a0  in  the  grazing  region.  Then,  analysis  on  the  only  available 
"shadow"  region  data  will  be  discussed.  Finally,  a  review  of  previous  terrain  models  and 
why  they  were  rejected  in  favor  of  a  sea  model  will  be  done. 


2.2  Calculation  of  Backscatter  Coefficient  (<j°) 

Because  many  factors  affect  the  magnitude  of  the  backscatter,  it  is  best  to  start  with 
the  basic  radar  range  equation  and  work  to  the  backscatter  measurement.  While  deriving 
the  equations,  it  will  become  clear  why  a  probabilistic  generation  of  g°  would  be  better 
than  a  deterministic  model  calculation. 

Starting  with  an  antenna  radiating  in  free  space  at  a  wavelength,  X,  and  signal 
power,  Pt  transmitting  isotropocially,  the  power  density,  P,  or  the  power  per  unit  area,  at  a 
range  R  from  the  antenna  is  calculated  to  be: 


P- 


P. 

AnR2 


(3) 


Since  radars  use  directive  antennas  to  increase  range,  antenna  gain,  Gt,  is  multiplied 
by  P  to  get  the  power  density  for  a  directed  antenna,  Pg.  Gain  from  the  antenna  is  defined 


G,  =  An 


(K/\ 

A2 


Pg=PG, 


(4) 

(5) 


as: 


with  Ae  being  the  aperture  or  receiver  part  of  an  antenna.  The  power  incident  at  a  point  at 
distance  R  from  the  transmitting  antenna,  and  on  its  boresight,  is  [4]: 


Pr=PS  A 


(6) 


p=pg2L 
r  (4 kJr1 


(7) 


The  radar  cross  section,  (RCS)  also  known  as  o,  of  a  target  may  be  defined  as  "the 
area  intercepting  that  amount  of  power  which,  when  scattered  isotropically,  produces  an 
echo  equal  to  that  received  from  the  object"  [12].  Simplified,  it  implies  a  measure  of  power 
that  would  be  returned  by  a  perfect  reflector  at  that  range.  The  power  received  from  a 
target  with  a  RCS  of  a,  Pa,  can  be  calculated  by  multiplying  the  antenna  power  density  by 
the  power  density  of  a  target  radiating  isotropically: 


P  a 

4nR2  4 kR2 

assuming  a  monostatic  radar  where  the  transmit  and  receive 
power  density  at  the  receive  antenna  is  simply  calculated  as: 


(8) 


antennas  are  co-located.  The 


=  PcA 


(9) 


This  results  in  the  basic  radar  range  equation  below: 


PtG2A2 
(4k)3  R4 


(10) 


where  Pr  =  Received  Power,  Pt  =  Transmitted  Power,  G  =  Antenna  Gain,  X  =  wavelength, 
R  =  range  from  antenna  to  target,  and  a  =  the  target's  radar  cross  section  (RCS)  at  that 
wavelength  [4].  Because  the  available  measured  data  is  overwhelmingly  monostatic,  the 
bistatic  case  is  ignored. 
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The  RCS,  a,  is  usually  normalized  for  an  ideal  target  or  surface  and  has  other 


dependencies  other  than  wavelength  [4].  This  will  discussed  in  section  2.3. 


Using  the  radar  range  equation  and  integrating  over  an  area  A  containing  N 
scatterers: 

Oj  i  =  1, 2,  3, . . .,  N  scatterers 

and  narrowing  down  the  beamwidth  such  that  it  approaches  0,  a  becomes  a0  and  the 
integral  collapses  to  a  single  system  constant  K  [1]. 

PN(e)=K(7°  (12) 

This  constant  can  be  calculated  for  each  testing  radar  system  to  allow  for  continuity 
between  separate  tests  and  radar  systems.  If  the  calculation  is  not  performed,  and  K  not 
compensated,  an  error  referred  to  as  Registration  error  is  introduced.  The  only  way  to 
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remove  this  error  is  to  have  a  complete  listing  of  all  of  the  measuring  radar’s  system 
parameters. 

Both  RADLAB  and  ARL  attempted  K  calculations  to  remove  errors  [13].  It  is 
assumed  that  all  registration  errors  between  the  RADLAB  data  and  ARL  data  equal  zero. 
For  future  data  measurements,  however,  at  least  a  factor  of  K  should  be  included  with  the 
published  data  for  confirmation  of  registration  error  removal. 

2.3  Terrain  Dependence  of  a0 

Since  the  constant  K  of  Equation  (12)  can  be  calculated  and  fixed,  the  return  power 
depends  on  c°.  The  amount  of  variance  of  c°  for  most  terrain  types  makes  it  very  difficult 
to  create  a  single  deterministic  equation  and  ultimately  is  the  reason  behind  the 
probabilistic  modeling  approach  that  is  used. 

The  value  of  o°  is  dependent  on  surface  moisture  content  as  demonstrated  by 
RADLAB  [7].  The  higher  the  moisture  that  is  in  the  surface,  the  higher  the  measured  c°. 
With  a  change  in  soil  moisture  from  desert  conditions  to  very  moist  ground,  RADLAB 
measured  a  a0  change  of  5-7  dB  [7]. 

In  the  shadow  region,  however,  the  a°  dependence  on  the  incidence  angle  becomes 
the  dominating  factor.  The  change  in  c°  varies  between  5-20  dB  from  70°  to  90°.  Below 
70°,  however,  o°  flattens  out  to  a  small  change  in  dB  until  0  =  30°  [1,3]. 

The  interaction  between  a0  (or  the  shape  of  the  overall  a0  curve)  appears  to  be 
dependent  on  the  surface  roughness.  If  the  surface  is  mostly  smooth  (snow,  asphalt,  dirt), 
the  curve  plateaus  and  then  quickly  drops  as  the  incidence  angle  is  increased  past  70°.  If 
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the  surface  is  rough  (trees,  grass,  rocks),  then  the  curve  will  gradually  drop  as  the  incidence 
angle  is  increased.  In  one  terrain  case,  urban  areas,  the  curve  remains  constant,  or  drops  a 
little,  no  matter  what  the  incidence  angle  is.  This  is  expected  due  to  the  number  of  comer 
and  angle  reflectors  that  can  be  found  in  buildings  and  structures  [2]. 

2.4  Simulated  Surface  Model 

As  stated  above,  a  terrain’s  roughness  appears  to  have  a  large  affect  on  the  shape  or 
percentage  change  in  o°  at  low-grazing  angles.  Meanwhile,  the  terrain  composition 
appears  to  affect  only  the  overall  value  of  a0  at  a  given  incidence  angle  (or  "height"),  not  its 
relationship  to  its  neighboring  values  (or  "shape  of  the  curve").  By  properly  modeling 
these  effects,  an  innovative  method  is  developed  for  modeling  G°  at  all  incidence  angles. 

If  the  curve  has  the  same  percentage  drop  in  values  for  terrain  of  the  same 
roughness,  but  not  of  the  same  surface  composition  or  material  makeup,  then  a  model  for 
one  surface  type  could  be  substituted  for  another  surface  type.  This  theory  can  tested  by 
creating  random  surfaces  of  varying  heights  and  assigning  electrical  constants  (permittivity 
and  conductivity)  to  the  surface  facets  [14]. 
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0  0 

Figure  4.  Random  Surface 


2.4.1  Rough  Surface. 

For  rough  terrain,  as  the  incidence  angle  changes,  the  change  in  a0  is  primarily 
dependent  on  two  effects,  geometric  changes  and  reflectivity  changes  [13].  Geometric 
changes  can  be  measured  in  the  number  of  facets  that  are  visible  for  a  given  incidence 
angle  and  the  average  slope  of  the  facets.  Reflectivity  changes  can  be  modeled  by  Fresnel 
reflection  coefficients  calculated  for  each  facet  as  given  by  Beckmann  [14]. 

Each  visible  facet  can  be  compared  to  a  scattering  point.  As  the  number  of  visible 
facets  increases  with  each  change  in  incidence  angle,  the  number  of  scatterers  increase, 
increasing  overall  c°. 

The  average  slope  of  the  facets  can  be  compared  to  the  amount  of  energy  reflected 
back  to  the  receiver.  As  facet  face  approaches  to  normal  incidence  with  the  receiver,  the 
amount  of  energy  reflected  back  to  the  receiver  increases.  Likewise,  as  the  facet 
approaches  grazing,  the  energy  reflected  back  decreases,  reducing  o°.  The  Fresnel 
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reflection  coefficients  calculate  the  amount  of  energy  each  facet  receives  and  with  the 
average  slope,  can  calculate  the  amount  of  energy  reflected  back  to  the  receiver  [14]. 


Figure  5.  Scattering  off  of  Random  Surface  [14] 


2.4.2  Fresnel  Reflectivity  Coefficients. 

The  coefficients  depend  on  intersect  angle,  the  electrical  properties  of  the  surface, 
and  the  polarization  of  the  intersecting  wave  [14]. 

n+  y2cost?-VT2 -sin2#  1n]  . 

R  = - .  —  Vertical  Polarization  (13) 

y2costf+VT2  -sin2$ 

R-  =  Horizontal  Polarization  (14) 

costf+VT2  -sin2i) 


y  =  —  Normalize  admittance 

\Vre 


(15) 


The  geometric  and  reflectivity  effects  on  a0  can  be  modeled  in  the  following  way: 
First,  create  a  random  rough  surface  with  a  given  average  height.  Then,  for  incidence 


angles  from  60°  <  0  <  90°,  calculate  the  number  of  facets  that  are  visible,  the  average  slope 
of  each  facet,  and  the  vertical  and  horizontal  Fresnel  reflectivity  coefficients  for  each  facet. 
Finally,  calculate  a  total  “a0”  from  the  coefficients,  number  of  facets,  and  average  slope. 


2.4.3  Random  Surface. 

For  this  set  of  graphs,  the  normalized  admittance,  Y,  is  calculated  for  grass, 


Y  =  2.24  [4], 


The  results  are  as  follows: 


Grass,  facets  visible  dB 


Figure  6.  Visible  Facets  in  dB  (Grass) 
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Coefficient  Slope  Angle  in  degrees 


60 


Grass,  Visible  Facet  Average  Slope  Angle 


Figure  7.  Average  Visible  Facet  Slope  Angle 
(Grass) 


Grass,  Visible  Facet  Vertical  Reflection  Coefficient 


Figure  8.  Vertical  Coefficient  (Grass) 
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Grass,  Visible  Facet  Horizontal  Reflection  Coefficient 


The  overall  values  of  each  graph  are  not  of  concern  in  this  study.  The  interaction  of 
the  values  for  each  surface  determines  the  value  of  “a0”.  By  calculating  the  percentage 
change  of  each  factor  (visible  facets,  average  slope  angle,  &  Fresnel  coefficients),  a  pattern 
emerges  that  shows  which  factor  changes  the  most  and  thus  has  the  greatest  effect  on  “c°” 
in  the  shadow  region. 

A  graph  of  a  theorized  “a0”,  calculated  for  each  angle  as: 
a0  =  (#  of  visible  facets)*(Average  facet  slope)  *(Vertical  Coeff.  +  Horizontal  Coeff.)  (16) 
allows  for  this  analysis. 
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Grass,  Total  Coefficent  dB 


It  is  interesting  to  note  that  the  curve  shapes  match  very  closely  to  the  “a0”  curve 
shapes  in  nature.  For  a  relatively  smooth  surface  (Low),  the  curve  drops  quickly.  This 
curve  shape  has  been  seen  in  measurements  of  roads  and  flat  ground.  The  rougher  the 
surface,  the  flatter  and  less  steep  the  curve  becomes,  as  seen  in  the  three  other  surfaces. 
This  curve  has  been  noticed  in  measurements  of  rougher  surfaces,  rocky  ground,  rough 
seas,  and  trees. 

The  effect  of  each  factor  is  measured  and  compared  to  the  other  factors  by  taking  a 
percent  difference.  Each  percentage  is  calculated  by  taking  the  difference  between  the 
lowest  and  highest  simulated  measurements  and  dividing  by  a  fixed  value.  The  fixed  value 
is  the  largest  value  each  simulated  measurement  could  obtain.  For  facets,  at  most  10,000 
possible  facets  could  be  seen.  For  average  slope,  the  value  is  90°.  For  the  Fresnel 
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Coefficients,  the  value  when  the  surface  is  normal  to  the  radar  beam  is  max.  a0  total  is  the 
exception  with  a  simple  percentage  change  of  (max-min)/max.  The  results  follow  in 
Table  1. 


Low 

Low-Medium 

Medium 

High 

A%  Visible  Facets 

79.9  % 

35.7  % 

22.1  % 

11.7% 

2.8% 

4.3% 

4.5% 

3.4% 

A%  Vert.  Coeff. 
max  mag.  (0.39) 

38.0  % 

15.0  % 

9.1  % 

5.0% 

A%  Hori.  Coeff. 
max  mag.  (-0.39) 

19.4  % 

11.3% 

7.8% 

4.7% 

A%  a0  Total 

83.5  % 

41.2% 

26.1  % 

14.9  % 

Table  1.  A%  Change  in  factors  (Grass) 


For  each  of  the  surface  heights,  the  percent  change  in  visible  facets  is  the  greatest 
and  appears  to  be  the  primary  effect  on  the  percent  change  in  total  c°.  If  this  pattern 
continues  for  other  values  of  Y  as  well  as  for  different  type  of  surface  shapes  (sea  surfaces), 
it  would  imply  that  the  theory  that  a0  is  mostly  affected  by  surface  roughness  at  low- 
grazing  angles  would  be  correct. 


2.4.4  Random  Shrubs  &  Trees. 

The  next  set  of  graphs  is  for  Y  calculated  for  shrubs  (small  average  random  height) 
and  trees  (large  average  random  height),  Y  =  3.16.  Again,  the  results  are  similar  as  before. 


Slope  Angle  in  degrees  Visible  Facets  dB 


Trees,  facets  visible  dB 


Figure  11.  Visible  Facets  in  dB  (Trees) 


Trees,  Visible  Facet  Average  Slope  Angle 


Figure  12.  Average  Visible  Facet  Slope  Angle 
(Trees) 
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Coefficient  Coefficient 


Trees,  Total  Coefficent  dB 


Figure  15.  Total  Coefficient  in  dB  (Trees) 


The  results  in  Table  2  are  calculated  the  same  way  as  Table  1 . 


Low 

Medium 

High 

A%  Visible  Facets 

77.5  % 

38.9  % 

25.1  % 

12.1  % 

A%  Ave.  Slope 

2.8% 

4.8% 

4.9% 

3.5% 

A%  Vert.  Coeff. 
max  mag.  (0.52) 

37.7  % 

14.4  % 

9.3% 

5.1  % 

A%  Hori.  Coeff. 
max  mag  (-0.52) 

19.3  % 

10.9  % 

7.9% 

4.8% 

A%  Total 

89.1  % 

43.5  % 

29.4  % 

16.3  % 

Table  2.  A%  Change  in  factors  (Trees) 


As  before,  in  the  simulated  grass  surface,  the  simulated  tree  surface  shows  same 
larger  percentage  change  in  the  number  of  visible  facets  compared  to  the  change  in  average 
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slope  and  the  coefficients.  This  supports  the  theory  that  “a0”  is  more  affected  by  surface 
roughness. 

2.4.5  Simulated  Sea  Surface. 

The  next  test  of  this  theory  is  to  see  if  a  non-random,  periodic  surface,  such  as  sea 
waves,  would  exhibit  the  same  behavior.  If  a  simulated  sea  surface  would  exhibit  similar 
changes  in  effects,  that  would  allow  use  of  a  sea  clutter  generator  to  model  the  overall 
shape  of  the  g°  curve  for  terrain  (ie.  the  change  in  a  from  one  angle  to  the  next). 

0.2 

0.15 

0.1 

0.05 

0 

60 

Figure  16.  Sea  Waves 

The  same  test  that  was  performed  on  the  random  surfaces  was  applied  to  the 
surface  shown  in  Figure  18.  The  results  are  as  follows: 
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incidence  angle  6 

Figure  18.  Average  Visible  Facet  Slope  Angle 
(Sea) 


Sea,  Visible  Facet  Vertical  Reflection  Coefficient 


Sea,  Total  Coefficent  dB 


Figure  21.  Total  Coefficient  in  dB  (Sea) 

The  sudden  spikes  in  the  graphs  of  facets  visible  and  total  “a0”  are  due  to  the 
surface  model  shape.  On  each  wave,  there  are  lines  of  facets  that  will  become  visible  at 
certain  viewing  angles.  Since  the  surface  has  ten  waves,  if  a  viewing  angle  is  passed,  10 
times  the  number  of  facets  will  become  visible  and  spike  the  line.  This  modeling  method 
does  not  affect  the  overall  data  or  its  conclusion. 


Low 

High 

A%  Visible  Facets 

62.3  % 

62.3  % 

6.2% 

6.2% 

1.1  % 

3.1  % 

0.9% 

1.2% 

A%  Vert.  Coeff. 
max  mag.  (0.79) 

45.2  % 

24.5  % 

13.7  % 

7.0% 

A%  Hori.  Coeff. 
max  mag.  (-0.79) 

10.6  % 

9.2% 

7.1  % 

5.0% 

A%  Total 

94.8  % 

65.1  % 

24.7  % 

15.2  % 

Table  3.  A%  Change  in  factors  (Sea) 
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Again,  since  the  %  change  is  largest  in  the  visible  facets,  for  all  four  surface 


heights,  this  shows  that  the  number  of  scatters  is  the  dominant  affect  on  a0. 

All  three  simulated  surfaces  (grass,  trees,  and  sea)  have  shown  the  same  effects  of 
number  of  facets  vs  the  average  slope  and  Fresnel  coefficients  on  “a0”.  This  implies  that  a 
model  that  works  on  sea  surface  “a0”  could  be  used  in  a  land  environment.  RADLAB 
noticed  a  similar  roughness  effect  and  that  will  be  discussed  next. 

2.4.6  RADLAB  Sensitivity  to  Roughness  and  Moisture  Content. 

RADLAB,  at  the  University  of  Michigan  modeled  the  effect  of  surface  roughness 

on  measured  data  and  compared  them  to  the  effect  made  by  moisture  content.  By  using  a 
soil  dielectric  model  given  in  Appendix  E  of  Ulaby  et  al.  [1],  they  measured  the  change  in 
G°  based  on  moisture  to  be  about  5.5  dB  at  10  GHz  and  less  than  3  dB  at  95  GHz.  This  was 
over  the  incidence  angle  range  55°  <  0  <  85°  [7].  This  implied  that  most  of  the  variation 
exhibited  by  0°  in  the  shadow  region  would  be  due  to  surface  roughness  since  c°  changed 
by  a  much  greater  dB  amount  over  the  angle  range. 

To  better  show  this  effect,  RADLAB  calculated  a  roughness  factor  for  6  surfaces 
they  had  data  at  95  GHz  [7].  Roughness  was  calculated  by: 

Roughness  =  ks  (17) 

where  k  =  2tc/X,  (wave  number)  and  s  is  the  rms  height  (Standard  Deviation).  The  smaller 
the  ks  is,  the  smoother  the  surface.  They  plotted  a0,  both  measured  data  and  theoretically 
calculated  curves,  for  two  types  of  surfaces,  smooth  surfaces  (ks  <  2)  and  rough  surfaces 
(ks  >  5)  [7].  The  next  two  figures  show  the  results. 
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coefficients  for  ks  <  2  [7] 


coefficients  for  ks  >  5  [7] 


For  smooth  surfaces  with  ks  <  2,  c°  decreases  at  a  steep  slope  up  to  88°  (the 
maximum  measured  data  angle).  In  contrast,  for  rough  surfaces  where  ks  >  5,  shows  a 
gentle  slope  relative  to  the  incidence  angle.  Both  graphs  show  the  need  for  a  curve 
generator  that  produces  shallow  slope  curves  for  very  rough  surfaces. 
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Because  the  three  previous  total  a  graphs  (random  grass,  random  trees,  and  water 
waves)  exhibit  the  proper  low-grazing  curves,  they  imply  that  surface  roughness  is  the 
primary  source  of  change  in  g°  in  the  shadow  region.  All  that  would  be  left  is  to  set  the 
overall  value  or  “height”  of  the  curves  such  that  the  generator  will  return  the  correct  a0. 
Gwm  GEN  can  provide  such  "height".  A  low-grazing  sea  model  for  use  in  curve  generation 
is  discussed  in  chapter  3. 

Now,  we  will  discuss  the  only  available  low-grazing  data  and  then  discuss  two 
other  low-grazing  terrain  models  that  were  considered,  but  finally  rejected. 

2.5  Analysis  of  Available  Grazing  Data 

The  only  available  data  for  grazing  angles  above  80°  was  from  the  Army  Research 
Laboratory  (ARL).  It  was  measured  in  support  of  the  Longbow  radar  program  from  1995- 
96.  They  used  a  pulse  radar,  measuring  at  95  GHz  (W-band),  with  an  effective  1° 
beamwidth.  The  incidence  angles  varied  from  85°  to  89°  with  most  of  the  data  measured 
near  85°  [6]. 

Because  only  the  mean  and  standard  deviation  of  the  measured  data  were  available, 
the  incidence  angle  is  assumed  to  be  85°.  Some  data  was  measured  closer  to  the  horizon 
and  could  induce  errors  in  the  published  data.  I  assumed  that  the  data  measured  closer  to 
the  horizon  did  not  deviate  from  the  85°  data  and  would  not  affect  the  published  mean  and 
standard  deviation  [6]. 

Another  possible  source  of  error  in  the  published  data  is  snow  on  the  test  surfaces. 
During  measurements  of  woods  and  plains-type  terrain  surfaces,  snowfall  covered  some  of 
the  test  areas.  While  the  snow  depth  was  measured  and  published,  the  areas  that  were 
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covered  while  measured  and  the  areas  that  weren't  covered  was  not  indicated  in  the  ARL 
report.  Since  the  number  of  snow  covered  measurements  was  minimal,  I  assumed  the 
published  data  was  not  greatly  affected  and  did  not  need  modification. 

A  number  of  surfaces  were  measured  by  the  Army.  However,  the  only  surfaces  that 
equate  to  available  data  from  RADLAB  were  grasses,  trees,  and  bushes  (short  vegetation). 
Some  of  the  trees  and  grasses  data  was  taken  with  snow  on  the  measured  cell.  As  stated 
before,  this  did  not  affect  the  use  of  the  published  data.  However,  it  must  be  noted  that 
RADLAB  has  separate  measurements  for  snow  [6]. 


Terrain  Cell  Type  (RADLAB  Designator) 

VV  Polarization 

Trees  (Coniferous)  ( B ) 

mean  (-17.6)  SD  (1.0) 

Trees  (Deciduous)  ( B  ) 

mean(-16.4)  SD(2.3) 

Grass  ( C ) 

mean(-18.3)  SD(1.7) 

Short  Vegetation  ( E ) 

mean(-19.3)  SD(1.9) 

Table  4.  ARL  measured  data,  Mean  a°  and 


Standard  Deviation  in  dB 


2.6  Currently  Available  Low-Grazing  Terrain  Backscatter  Models 

In  researching  backscatter,  two  land  clutter  models,  besides  the  squared  cos  term 
mentioned  in  Chapter  1,  were  repeatedly  discussed.  Georgia  Tech  and  RADLAB  at  UM 
developed  deterministic  models  to  fit  their  available  data.  While  neither  model  was  used  in 
the  final  solution  or  for  validation,  a  discussion  of  both  models  is  included  for  comparison 
to  the  final  developed  model. 

Both  models  suffer  from  two  common  problems.  The  first,  and  more  serious, 
problem  relates  to  their  deterministic  properties.  Both  models  return  a  mean  for  the  terrain 
constants  they  have.  However,  neither  model  generates  data  with  the  randomness  of 


30 


nature.  The  second  problem  depends  on  their  ability  to  accept  new  measured  data.  Since 
both  models  are  curve-fitting  algorithms,  any  new  measured  data  would  require  more 
curve-fitting  analysis  before  the  data  could  be  used 


2.6.1  Georgia  Tech  Land  Model. 

Because  of  the  difficulty  in  modeling  a  complex  terrain  surface,  personnel  at 
Georgia  Tech  created  an  empirical  model.  Instead  of  measuring  the  effects  of  surface 
roughness,  surface  composition  (moisture  content),  distribution  of  scatterers,  and  other 
factors,  the  model  was  fitted  to  a  selected  subset  of  available  data  [3]. 

-d/H 

A,  B,  C,  &  D  are  derived  constants  based  on  fitting  the  curve  to  the  measured  data, 
not  the  measured  frequency,  surface  type  and  roughness,  or  other  effects.  Gh  is  the  standard 
deviation  of  the  surface  in  decimeters  and  X  is  the  beam  wavelength. 

While  there  are  available  constants  for  almost  all  of  the  nine  terrain  surfaces,  the 


G°  =  A(0  +  Cf  exp 


Georgia  Tech  model  does  not  have  the  frequency  range  (only  in  Ka  and  W-band)  and  angle 
range  (35°  <  0  <  70°)  that  is  desired.  Also,  in  comparison  to  the  sea  curves,  the  land  model 
does  not  drop  near  the  horizon  as  usually  seen  in  data  (see  Figure  24). 


o°  Calculated,  GT  land  Model  vs  Sea  Model  for  W-band,  Trees 


Figure  24.  Georgia  Tech  Land  Model  vs.  Sea 
Clutter  Model,  GT  Curve  (Top),  Sea  Clutter 
Model  (Bottom) 


Since  the  Georgia  Tech  model  was  developed  for  "plateau"  region  operations 
(35°  <  0  <  70°),  its  shadow  region  operations  are  limited  at  best  (as  seen  in  Figure  24).  It  is 
for  this  reason  that  this  model  was  not  considered  for  final  usage  with  Owm  GEN  . 


2.6.2  RADLAB  Model. 

Personnel  at  the  University  of  Michigan  base  this  model  on  the  extensive  collection 
of  correlated  measured  data  [1].  Like  the  Georgia  Tech  model,  the  RADLAB  model  is 
fitted  to  a  selected  subset  of  measured  data.  The  RADLAB  model  contains  two  equations: 
a  mean  o°  equation  and  a  standard  deviation  equation  based  on  the  measured  data.  The 
idea  is  that  a  more  random  G°  could  be  generated  using  these  two  equations. 
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m  =  Pl+P2  exp(-  P3d)+  P4  cos {p5d  +  P6 ) 


(19) 


SD(Q)=  +  M 2  exp(- M30)  (20) 

where  Px  and  Mx  are  frequency  and  terrain  dependent  constants. 

While  this  model  and  Cwm  GEN  are  based  on  the  same  collection  of  data,  the 
RADLAB  model  is  more  limited  than  Owm  GEN  in  the  frequencies  and  terrain  types  that 
are  available.  This  limitation,  as  well  as  RADLAB's  inability  to  accept  new  data  without 
curve-fitting,  made  the  RADLAB  model  less  desirable  to  use. 

2.7  Summary 

The  overall  problem  of  measuring  and  simulating  c°  for  terrain  was  shown.  The 
measurements  of  the  faceted  random  and  sea  surfaces  appear  to  suggest  the  theory  that  at 
low-grazing  angles,  most  of  the  effect  on  c°  is  due  to  the  overall  surface  roughness,  not 
surface  composition.  This  allows  for  a  low-grazing  sea  model  to  be  use  to  model  terrain 
clutter.  The  ARL  data  was  analyzed  and  discussed  for  possible  errors.  It  was  decided  that 
the  potential  errors  would  have  very  little  effect  on  the  final  product  and  could  be  safely 
disregarded.  Finally,  two  currently  used  land  models  were  compared  and  rejected  due  to 
their  operational  limits. 
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III.  Methodology 


3.1  Introduction 

The  simulated  surface  example,  set  forth  in  section  2.4,  allows  for  the  usage  of  a  o° 
generator  of  a  surface  to  be  used  with  a  surface  with  a  similar  roughness.  All  that  would  be 
needed  to  be  determined  would  be  the  overall  value  of  the  curve.  Gwm  GEN  completes 
this  task. 

As  seen  in  the  last  chapter,  there  are  very  few  terrain  models  available,  and  of  those, 
none  have  the  frequency  range,  terrain  mixture,  and  angle  set  that  is  needed.  However, 
there  are  two  sea  clutter  G°  models,  Georgia  Tech  and  the  Center  of  Naval  Analysis  (CNA), 
which  contain  the  angle  set  and  dependence  on  wave  height  (surface  roughness)  that  is 
desired  [5].  Because  the  Georgia  Tech  model  had  the  largest  available  frequency  range,  it 
is  the  selected  model  and  used  in  the  final  generator. 

This  chapter  will  discuss  the  development  of  the  Georgia  Tech  model  and  describe 
how  it  is  integrated  with  Owm  GEN.  Then  the  completed  model  will  be  compared  to  Army 
data  in  four  cases  where  the  Army  data  coincides  with  the  RADLAB  data  (basis  for  Gwm 


GEN). 


3.2  Sea  Model  Theory 


The  model  theory  is  based  on  Currie’s  description  of  the  millimeter  portion 
(v  >  30  GHz)  [3].  The  original  model  was  developed  and  described  by  Horst  and  Perry, 
however  the  original  document  is  not  available  [15].  The  greater  than  millimeter  wave 
portion  (v  <  30  GHz)  was  derived  from  the  ALARM  91  model  from  SAIC  [5]. 


Symbol 

Definition 

0 

Incidence  angle 

X 

Radar  wavelength 

<t> 

angle  between  boresight  and  wind  direction 

V  wind 

Velocity  of  wind  (m/s) 

HHorVV 

Polarization  of  transmitted  and  received  waves 

Table  5.  Input  parameters  for  Sea  Clutter  Model 


Symbol 

Definition 

Variable  Dependence 

h\vave 

Wave  height 

Vwjnd 

r 

roughness  parameter 

K  (J),  hwave 

Ai 

multipath  interference  factor 

r 

Au 

upwind-downwind  factor 

0,  x,  ux 

Ui,  U2,  u3 

est.  parameters  based  on  meas.  data 

Aw 

Wind  speed  factor 

Vwind,  qw 

qw 

power  factor 

X,  qwi 

qwi 

est.  parameter  based  on  meas.  data 

_o 

<7  XX 

calculated  backscatter  for  HH  or  VV 

X,  0,  Aj,  Au,  Aw,  cx,  zx 

Ci,  C2,  Zi,  Z2,  Z3,  Z4,  Z5,  Z6 

est.  parameters  based  on  meas.  data 

Table  6.  Internal  variables  for  Sea  Clutter  Model 
and  dependence 


The  problem  of  developing  a  model  for  sea  backscatter  is  a  complicated  one  at  best, 
since  it  is  dependent  upon  radar  frequency,  sea  state,  incidence  angle,  wind/sea  direction. 


polarization,  and  other  factors.  Georgia  Tech  researchers  designed  an  analytical  model  to 
describe  sea  clutter  backscatter  using  empirical  constants  derived  from  measured  data  [3]. 

The  clutter  model  calculates  a0  as  the  product  of  three  variables:  multipath 
interference,  sea  direction,  and  wind  speed.  Each  of  these  factors  is,  in  turn,  a  function  of 
the  appropriate  variables. 

The  multipath  interference  effect  factor,  between  the  direct  and  scattered  fields 
from  the  surface  for  horizontal  polarization,  is  derived  from  forward  scatter  theory 
assuming  a  Gaussian  distribution  of  surface  height  with  a  standard  deviation  of  g°hh 
A  roughness  parameter,  r,  is  defined  to  be: 

r  =  (14.4A  +  5.5)^  (21) 

A 

where 


=  0.3481 


(  V  V 1 

1  V  wind  1 

11 


(22) 


The  average  wave  height  (hwaVe)  is  defined  to  match  the  sea-state  descriptors  as  set 
in  Skolnik  [9,  pi 36]. 


Sea  State 

Wind  Speed  (m/s) 

Wave  Height  (m) 

1  (smooth) 

<3.5 

<0.3 

2  (slight) 

3.5 -6.2 

0.3 -0.9 

3  (moderate) 

6.2 -8.2 

0.9 -1.5 

4  (rough) 

8.2 -9.8 

1.5 -2.4 

5  (very  rough) 

9.8-11.8 

2.4 -3.7 

6  (high) 

11.8-15.4 

3.7-6 

7  (very  high) 

15.4  -  23.2 

6-12.2 

Table  7.  Sea-State  Surface  Descriptors  (Numbers 
are  not  integer  values  due  to  wind  speed  set  in 
knots  and  wave  height  set  in  feet.) 
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For  this  model,  the  sea  states  do  not  need  to  be  integer  values.  A  sea  can  be  described  to  be 
3.5  if  the  wind  speed  and  wave  height  match  for  that  particular  sea  state.  In  the  final 
integrated  model,  a  number  of  surfaces  are  assigned  non-integer  sea-states  for  roughness 
parameters.  More  information  can  be  found  in  Chapter  4. 

The  interference  term  A*  is  given  by 


(23) 


The  behavior  of  a0  versus  incidence  angle  falls  into  two  distinct  regions:  a  low 
grazing  angle  region  in  which  G°  is  a  strong  function  of  angle,  and  a  plateau  region  in 
which  G°  is  approximately  independent  of  the  angle.  Due  to  this,  only  <T°  data  from  the 
plateau  region  was  used  in  determining  the  parametric  dependencies  for  the  model. 

The  sea  direction  term  is  based  on  upwind-downwind  data,  although  the  reference 
for  the  aspect  angle  should  be  the  sea  wave  propagation  vector  rather  than  the  wind 
vector.  The  up-down  data  is  dependent  on  incidence  angle  0,  increasing  as  0  approaches 
the  horizon,  but  approaching  a  finite  value  at  horizon  [3].  Since  there  is  insufficient  data  to 
determine  clearly  the  functional  form  of  the  dependence,  a  simple  function,  cos  <|),  is 
chosen,  where  <f>  is  the  angle  between  antenna  boresight  and  the  upwind  direction  [3].  For 
system  simplicity,  <|)  is  fixed  at  90°,  but  can  be  changed  if  necessary.  The  upwind- 
downwind  interference  term  Au  is  defined  to  be: 


\  =  exp[w,  cos  0(1  -  2.8 0)(A  +  w2)"3  ] 


(24) 


Currie  states  that  the  dependence  of  a0  on  the  sea  state  has  been  shown  to  be  more  a 
function  of  wind  speed  than  of  wave  height  [3].  This  is  interesting  because  the  wave  height 
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(hWave)  is  calculated  from  the  wind  speed.  For  simplicity,  the  wind  speed  term,  Aw,  is 
defined  as: 


Av  = 


(25) 


qw  =  qwx  (A  +  0.2)~° 4  (26) 

Through  curve  fitting  of  the  available  measured  data,  the  final  c°  equation  was 
found  to  be: 

o°HH  =  c.Afl'MAA  (V  <  10  GHz)  (27) 

and 

Kh  =  c^AAA  (v  >  10  GHz)  (28) 

The  correction  factor  for  vertical  polarization  has  been  found  to  be  a  linear  fit 
compared  to  the  c°  for  horizontal  polarization.  Using  the  data,  c° w  is: 

<  =°hh  -{hwave+zx)zHe+z,y*xnw*  (29) 

where  cx  and  zx  are  empirically  derived  frequency  constants. 

Cross  polarization  (a°hV  &  ct°Vh)  data  was  not  available  when  the  model  was 
developed.  Ulaby  et  al.  have  found  that  near  the  horizon,  measured  o° w,  o°hV,  &  0%  are 
within  2.0-4.0  dB  of  each  other  with  cross-polarization  usually  less  than  linear-polarization 
[7].  The  final  model  will  use  a°vv  curves  for  the  cross-polarizations. 

The  final  result  of  this  model  are  sea  clutter  curves  as  shown  in  the  next  figure.  As 
the  sea  surface  becomes  rougher,  the  curve  shape  flattens  and  does  not  drop  off  until  very 
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close  to  the  horizon.  This  effect  has  been  shown  in  measured  data  [4,  7]  and  is  the  desired 


behavior  for  the  final  integrated  model. 


Sea  Model  Curves  X-band 


Figure  25  Sea  Model  Curves  for  5  Different  Sea 
States 


In  the  final  integrated  model,  an  "elbow"  effect  will  occur  shortly  after  the 
integration  point  between  the  measured  data  and  the  sea  curves.  This  is  due  to  the  final 
dropping  off  of  the  sea  model.  Validation  tests  in  section  3.5  will  show  that  the  "elbow" 
effect  has  no  adverse  consequences  on  the  final  model  results. 
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3.3  Wilson/Mediavilla  Model 


The  following  equation  is  the  basis  behind  Gwm  GEN: 

<7  =/?•(- In  (£/)“)  + 7  (30) 

The  oc-p-y  combinations  are  set  such  that  they  would  yield  the  required  mean  g°  and 
SD(0)  values  given  by  the  RADLAB  data.  A  typical  set  of  values  would  be  of  the  form: 


a  =  0.06, 0.061,0.062, ... ,  2.0 

(31a) 

P  =  0.01, 0.01 1,  0.012,  0.013, ... ,  14.0 

(31b) 

y=0.01 

(31c) 

When  a  combination  of  a,  p,  and  y  generates  the  mean  G°  and  SD(0)  for  a  given 
surface  and  angle,  those  values  of  a,  P,  and  y  are  placed  in  a  lookup  table.  When  a  G°  is 
requested,  the  model  retrieves  the  necessary  a,  p,  and  y  and  combines  them  with  a  U(0,1) 
randomly  generated  number.  The  method  of  generating  a,  P,  and  y  for  a  given  terrain  and 
angle  is  described  in  [2]  and  the  reader  is  invited  that  work  for  more  information.  For  this 
system,  the  model  is  representative  of  real  terrain  scattering  G°  measurements  for 
0°  <  0  <  80°. 

3.4  Model  Integration 

Since  Owm  GEN  is  based  on  measured  data,  the  integrated  model  will  return  the 
Gwm  GEN  g°  whenever  possible.  The  sea  model  is  only  implemented  when  the  requested 
angle  is  unavailable  through  Gwm  GEN. 
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If  such  an  angle  is  requested,  the  system  inquires  Gwm  GEN  for  the  highest 
measured  data  incidence  angle.  g0,s  values  are  calculated  for  the  highest  angle,  and  1°  and 
2°  below  the  highest  angle.  These  values  are  averaged  to  create  an  artificial  “plateau”  for 
the  sea  model  to  start  calculations.  This  plateau  is  usually  calculated  around  0  =  75°  for 
most  terrain  types  and  frequency  bands. 

This  method  of  plateau  calculation  combines  the  probabilistic  nature  of  Gwm  GEN 
with  the  accuracy  and  stability  that  comes  from  averaging  a  set  of  points  [2,  13].  All 
previous  near-horizon  models  provide  the  same  answer  for  a  given  terrain  type  and 
frequency  band.  By  using  Gwm  GEN,  the  probabilistic  nature  of  is  harnessed  while  the  sea 
model  gains  a  stable  point  to  start  from. 
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The  probabilistic  properties  of  Gwm  GEN  pass  though  the  sea  model  to  affect  any  low- 
grazing  angle  calculated  G°.  This  property  is  shown  in  the  four  validation  figures  in  the 
next  section. 

3.5  Model  Validation 

Since  the  Army  data  was  measured  in  W-band,  only  four  terrain  types  match  to 
RADLAB  data.  These  terrain  types  are  grasses,  vertical-vertical;  coniferous  trees 
(evergreens),  vertical-vertical;  deciduous  trees  (with  leaves),  vertical-vertical;  and  short 
vegetation,  vertical-vertical. 

The  single  vertical  line  at  85°  is  the  range  of  the  Army  data  out  to  one  standard 
deviation.  The  two  lines  running  outside  of  the  centerline  (the  mean  o°)  are  the  minimum 
and  maximum  o°’s  calculated  by  the  final  integrated  model. 

As  mentioned  before  in  section  3.2,  each  figure  shows  an  "elbow"  effect  at  some 
point  after  the  merger  point  from  measured  data  to  sea  curve  (around  0  =  80°).  That  is,  the 
resulting  curve  drops  faster  than  it  did  at  lower  angles.  The  following  figures  show  this 
effect  resulting  in  the  final  model  matching  the  measured  data.  This  is  most  noticeable  in 
Figure  33  (Short  Vegetation)  where  the  ARL  measured  data  line  (0  =  85°)  is  lower  than  the 
RADLAB  measured  data  (0  <  80°). 
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With  the  exception  of  short  vegetation,  the  final  model  c°  mean,  maximum,  and 
minimum  are  within,  or  very  near,  the  first  moment  of  the  measured  data.  The  disconnect 
of  min/max  of  the  final  model  for  short  vegetation  to  the  first  moment  is  due  to  the 
apparent  small  sample  size  of  the  ARL  data.  However,  the  mean  of  the  final  model 
matches  closely  to  the  mean  of  the  measured  data. 

3.6  Summary 

This  chapter  has  discussed  the  operation  of  the  sea  model  used  and  how  it  interfaces 
with  CT\vm  GEN  to  create  the  final  system  model.  A  brief  discussion  of  Owm  GEN 
followed  and  finally,  validation  of  the  final  model  was  done  with  measured  data  from  ARL. 
Chapter  4  has  test  runs  of  the  final  integrated  model  over  all  nine  terrain  surfaces. 


47 


IV.  Results  and  Discussion 


4.1  Overview 

Due  to  the  lack  of  measured  data  outside  of  the  four  cases  used  to  validate  the  final 
model,  any  discussion  is  limited  to  how  the  model  performs  and  if  the  different  polarization 
cases  (HH,  W,  and  HV)  relate  to  each  other  as  they  do  in  measured  scenarios. 

The  following  chapter  displays  model  results  for  each  terrain  type  in  two  different 
frequency  bands.  Comments  are  made  for  each  test  run  as  to  the  wave  height  used  and  if 
the  polarizations  behave  as  theorized. 

4.2  Sea  State  decisions 

The  basis  for  the  decision  of  each  sea-state  assignment  to  each  terrain  type  is  very 
simple.  Whenever  measured  data  near  the  horizon  was  available,  the  sea  state  was  chosen 
experimentally  by  trial  and  error.  Since  only  three  RADLAB  terrain  types  (grasses,  trees, 
and  short  vegetation)  were  available  under  the  ARL  data,  only  three  sea  states  were 
determined  this  way.  The  other  six  terrain  type  sea  states  were  determined  by  comparing 
the  three  experimentally  determined  terrain  types  with  the  other  terrain  types.  Surfaces 
smoother  than  grass  (roads  and  snow)  were  assigned  smaller  sea  states.  Surfaces  rougher 
than  trees  (Urban  areas)  were  assigned  larger  sea  states.  All  other  surfaces  were  assigned 
sea  states  in  between. 

After  all  terrain  types  were  assigned  sea  states,  initial  graphs  were  created  to  check 
if  the  measured  data  had  curvature  that  did  not  match  with  the  initial  sea  state.  If  that 
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happened,  the  sea  state  was  either  slightly  raised  or  lowered  until  the  sea  curve  matched 
with  the  measured  data  curve. 

4.3  Model  Results  for  selected  tests 

Before  discussing  the  final  model  results,  the  behavior  of  the  polarizations  needs  to 
be  recapped.  Measurements  performed  by  RADLAB  in  W-band  grazing  angle  regions 
have  shown  that  VV  and  HH  polarization  measurements  were  within  2-5  dB  of  each  other. 
Meanwhile,  HV  polarization  measurements  were  usually  10-15  dB  below  VV  and  HH  with 
VH  measurements  equaling  HV  [7]. 

In  the  figures  below,  the  model  data  fits  the  behavior  that  RADLAB  noted  in  most 
cases.  In  those  figures  where  the  polarizations  did  not  match  RADLAB  observations,  the 
measured  data  sample  size  was  noticeably  low.  Mediavilla  noted  that  measurement  cases 
with  small  sample  sizes  would  result  in  erroneous  results  more  often  than  surfaces  with 
high  number  of  sample  measurements  [2].  All  surfaces  with  very  low  sample  sizes  and 
very  large  erroneous  results  are  identified  by  Mediavilla  and  are  not  used  in  the  final 
integrated  system  model. 
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4.2.1  Dirt  &  Rocks  (RADLAB  Surface  A). 

The  sea  state  chosen  for  this  was  0.8.  This  was  chosen  based  on  the  grass  sea  state 

of  1.0.  Since  a  soil  and  rock  field  is  just  smoother  than  a  field  of  grass,  this  decision  is 

justified. 

In  the  final  results,  the  L-band  curves  are  smooth  with  the  HH  and  VV  polarizations 
within  a  few  dB  and  the  HV  polarization  about  10  dB  lower.  The  Ka-band  curves  show  a 
quick  dropping  off  at  80°.  The  results,  however  are  still  valid. 


o°  Calculated,  L-band,  Dirt  &  Rocks 
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Figure  34.  L-Band 


50 


a"  Calculated,  Ka-band,  Dirt  &  Rocks 
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Figure  35.  Ka-band 
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4.2.2  Trees  (RADLAB  Surface  B). 

The  sea  state  of  3.0  is  chosen  for  trees.  This  is  an  average  of  the  validation  results 
of  the  ARL  data  for  deciduous  and  coniferous  trees  (sea  states  4.0  and  2.0  respectifully). 

For  both  graphs,  the  polarizations  and  curvature  match  very  well  to  expectations  of 
HH  and  VV  above  HV  as  well  as  a  smooth  integation  and  drop  off.  The  linear  and  cross 
polarizations  are  not  seprated  by  10  dB,  but  this  is  thought  to  be  due  to  the  large  number  of 
vertical  and  horizontal  reflectors  (branches)  present  in  trees. 
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4.2.3  Grasses  (RAD LAB  Surface  C). 

Due  to  ARL  data  validation,  the  sea  state  of  grass  is  set  to  1.0.  This  does  not  mean 
that  all  other  sea  states  are  based  off  of  grass;  it  just  happened  to  work  out. 

The  curvature  for  both  graphs  are  mostly  smooth  with  the  exception  of  the  linear 
polarizations  in  Figure  39.  The  quick  drop  is  acceptable  and  not  seen  to  be  a  problem. 


cf  Calculated,  C-band,  Grasses 
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Figure  38.  C-Band 
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o°  Calculated,  Ku-band,  Grasses 
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Figure  39.  Ku-Band 
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4.2.4  Shrubs  (RAD LAB  Surface  D). 

Shrubs  sea  state  is  set  between  short  vegetation  and  trees  at  1 .7.  This  is  because 
shrubs  are  sometimes  considered  short  trees.  The  curvature  is  not  as  smooth  as  desired, 
however,  it  is  acceptable.  The  polarizations  behave  as  desired  with  cross  ploaziation  about 
7-10  dB  below  linear  polarizations. 
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4.2.5  Short  Vegetation  (RADLAB  Surface  E). 

Short  Vegetation  is  the  last  of  the  ARL  data  validated  sea  states.  This  sea  state  is 

set  at  1.3.  The  curvature  looks  good  with  a  small  drop  in  Figure  43.  Again,  this  drop  is 

acceptable.  The  polarizations  are  as  expected  with  a  10  dB  difference  in  L-band  between 

cross  and  linear. 

o°  Calculated,  Ku-band,  Short  Veg. 
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Figure  42.  Ku-Band 
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o°  Calculated,  L-band,  Short  Veg. 
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Figure  43.  L-Band 
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42.6  Roads  (RADLAB  Surface  F). 

Roads  sea  state  is  set  very  low  at  0.5  since  the  surface  is  usually  smooth.  The 
curvature  shape  matches  well  to  the  measured  data.  The  HH  polarization  is  lower  than 
expected,  (closer  to  HV  than  to  VV)  however  it  is  still  acceptable. 
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o°  Calculated,  Ka-band,  Roads 


4.2.7  Urban  Areas  (RADLAB  Surface  G). 

Urban  areas  are  only  included  for  experimental  purposes.  Mediavilla  feels  the 

measurements  sample  size  is  too  small  and  too  erratic  to  be  accurate  [2].  He  recommends  a 

larger  data  set  needs  to  be  measured  before  the  model  is  accuate.  The  sea  state  is  set  at  6.0 

but  can  be  changed  with  new  validation  data. 

□"  Calculated,  L-band,  Urban  Areas 
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Figure  46.  L-Band 
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c f  Calculated,  X-band,  Urban  Areas 
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Figure  47.  X-Band 


4.2.8  Dry  Snow  (RAD LAB  Surface  H). 

For  both  dry  and  wet  snow,  the  curvature  and  polarizations  perform  as  expected. 
The  sea  state  for  dry  snow  is  set  to  0.5  due  to  its  general  surface  smoothness. 
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4.2.9  Wet  Snow  ( RADLAB  Surface  I). 

As  stated  before,  the  curvature  and  polarizations  are  as  expected.  The  sea  state  for 
wet  snow  is  0.6,  slightly  higher  than  dry  snow.  This  is  due  to  the  higher  mositure  content 
of  wet  snow  which  has  been  shown  to  result  in  a  higher  a0  when  compared  to  dry  snow 


[1,7]. 


cf  Calculated,  S-band,  Wet  Snow 
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Figure  50.  S-Band 
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Angle  of  Incidence  6 

Figure  51.  Ka-Band 


V.  Conclusions  and  Recommendations 


A  complete  incidence  angle  G°  model,  Gwmoc  GEN,  for  use  in  terrain  modeling 
simulations  was  developed  and  tested. 

Having  measured  data  for  only  four  surface  types  at  one  frequency  makes  it  very 
difficult  to  validate  the  final  model.  Still,  this  is  the  first  a0  model  for  incidence  angles  0° 
<  0  <  90°  for  multiple  terrain  types  and  frequency  bands.  Any  new,  measured  data  will 
serve  only  to  improve  this  model. 

Military  applications  for  this  final  model  are  numerous.  This  model  was  originally 
developed  for  backscatter  removal  from  simulated  radar  returns  for  improved  target 
detection.  This  model  is  not  limited  to  backscatter  removal  but  can  also  be  used  if  the 
backscatter  is  the  target.  When  matched  with  DTED  data,  terrain  identification  can  assist 
with  ground  operations.  Finally,  surveillance  on  crops  and  forest  conditions  can  be  quickly 
completed,  giving  potential  dry  conditions  where  wildfires  could  become  a  threat.  Other 
uses  will  be  identified  as  the  model  improves. 

The  final  model  works  as  specified.  However,  there  are  several  ways  to  improve  it. 
Recommendations  for  future  works  include: 

1.  Perform  a  complete  remeasurement  of  the  nine  terrain  types  in  the  full  angular 
range  (0°  -  90°)  in  the  7  radar  bands  using  at  least  two  different  radars  in  the 
same  test.  The  use  of  two  different  radars  will  allow  for  the  removal  of  any 
constants  introduced  due  to  registration  errors. 
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2.  Introduce  a  propagation  model  to  the  low-grazing  portion  of  the  model.  A 
system  from  John  Hopkins  is  used  by  the  Navy  to  provide  site-specific  G° 
measurements  based  on  DTED  data  and  a  propagation  model.  The  system 
works  by  sending  a  series  of  pulses  at  an  area  and  processing  the  returns 
through  a  propagation  model  to  remove  any  atmospheric  effects.  Then,  using 
DTED  data  for  the  area  observed,  the  system  measures  the  <7°  returned  from  that 
terrain.  The  measured  a0  is  then  used  to  remove  backscatter  energy  from 
subsequent  pulse  returns.  This  system  suffers  from  problems  with  inaccurate 
weather  forecasts,  inaccurate  DTED  data,  and  the  possibility  of  extra  targets  not 
on  the  DTED  data  being  measured  and  excluded  [11].  A  compilation  of  this 
final  model  and  the  John  Hopkins  system  could  result  in  a  more  complete 
model. 

3.  Research  the  vector  cross-polarization  terms.  This  research  assumed  that  the 
cross-polarization  terms,  a°VH  and  a°Hv,  were  simple  polarizations.  However, 
there  are  vector  cross-polarization  terms  which  affect  how  the  transmitted  and 
received  energies  relate  [4, 13]. 
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Appendix  A.  Terrain  Backscatter  Code  in  MATLAB® 


A.l  Introduction 

The  following  two  functions  are  included  as  a  reference  for  the  figures  and 
equations  described  in  this  work.  The  WMOC  function  (Wilson,  Mediavilla,  and 
O’Connor)  is  the  main  function  that  calls  Wilson/Mediavilla  Weibull  function  and  the 
Seamodel  function.  The  reader  is  recommended  to  read  Mediavilla  [2]  on  the  creation  and 
format  of  the  Weibull  coefficient  table  that  is  the  basis  of  the  function. 

MATLAB®  was  chosen  as  the  programming  language  due  to  its  ease  and 
availability.  A  FORTRAN®  program,  based  on  these  functions,  is  available  from  the 
author  and  AFRL/SNAS. 


A.2  WMOC  Function 

function  sigmadb  =  wmoc(angle, frequency, terrain, polarization) 

%  WMOC  Wilson,  Mediavilla,  O'Connor  Terrain  Backscatter  Generator 
%  Written  23  Mar  99 

%  WMOC  will  return  a  backscatter  sigma  value  for  a  given  incidence 
%  angle  (degrees),  frequency  (1EEE/ANSI  STD  527),  terrain  type 
%  (RADLAB  nomenclature),  and  polarization  (HH,  VV,  HV) 

%  S  =  WMOC(angle,  frequency,  terrain,  polarization) 

%  Load  terrain  file  to  determine  highest  measurement  incidence  angle 
%  Nomenclature  e  (terrain  type)  w  (frequency  band)  w  (polarization) 

%  See  Mediavilla  for  terrain  file  creation  and  format, 
load  ewvv 

%  Determine  maximum  measured  incidence  angle 
max  angle  =  max(uldo(:,6)); 


if  angle  <=  maxangle 
%If  the  requested  angle  has  measured  data 
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%Generate  a  random  number 
U  =  rand(l); 


%Retrieve  alpha,  beta,  and  gamma  of  Weibull  distribution  for 
%the  requested  angle 
row  =  find(uldo(:,6)==angle); 
alpha  =  uldo(row,12); 
beta  =  uldo(row,13); 
gama  =  uldo(row,14); 

%Calculate  Sigma 

sigma  =  beta*((-log(U)).Aalpha)+gama; 

%Convert  to  dB  and  return 
sigmadb  =  20*logl0(sigma); 


else 

%If  the  requested  angle  does  not  have  measured  data 

switch  upper(frequency) 
case'L',  freq=1.5e9; 
case  'S',  freq  =  3.0e9; 
case  'C',  freq  =  6.0e9; 
case  'X',  freq  =  10.0e9; 
case  'KU',  freq  =  15e9; 
case  'KA',  freq  =  35e9; 
case'W',  freq  =  95e9; 

otherwise,  disp('Unknown  Frequency  Band.') 
end 

switch  upper(terrain) 

case  'A',  sea  =  0.8;  %  Soils  &  Rock  Surfaces 

case  'B',  sea  =  3.0;  %  Trees 

case  'C',  sea  =  1.0;  %  Grass 

case  'D',  sea  =  1.7;  %  Shrubs 

case  'E',  sea  =  1.3;  %  Short  Vegetation 

case  'F',  sea  =  0.5;  %  Roads 

case  'G',  sea  =  6.0;  %  Urban  Area 

case  'H',  sea  =  0.5;  %  Dry  Snow 

case  T,  sea  =  0.6;  %  Wet  Snow 

otherwise,  disp('Unknown  Terrain  Type.') 

end 

switch  upper(polarization) 
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case  'VV,  pole  =  0;  %  Vertical-Vertical 
case  'HH',  pole  =  1;  %  Horizontal-Horizontal 
case  'HV',  pole  =  0;  %  Horizontal-Vertical 
case  'VH',  pole  =  1;  %  Vertical-Horizontal 
otherwise,  disp('Unknown  Polarization  Type.') 
end 

%  Sea  model  does  not  work  at  horizon.  Value  at  horizon 
%  is  set  at  0.1  degrees  below  horizon 
if  angle  ==  90 

sigmadb  =  seamodel(89.9,  sea,  freq,  pole); 
else 

sigmadb  =  seamodel(angle,sea, freq, pole); 
end 

%  Set  measured  data  plateau  for  sea  model  to  start  from 

%  Calculate  first  point  of  measured  data  vs  sea  model 
%  Generate  a  random  number 
U  =  rand(l); 

%  Retrieve  alpha,  beta,  and  gamma  of  Weibull  distribution  for 
%  the  requested  angle 
row  =  find(uldo(:,6)==maxangle); 
alpha  =  uldo(row,12); 
beta  =  uldo(row,13); 
gama  =  uldo(row,14); 

%  Calculate  Sigma 

sigma  =  beta*((-log(U)).Aalpha)+gama; 

plateau(l)  =  20*logl0(sigma)  -  seamodel(maxangle,sea, freq, pole); 

%  Calculate  second  point  of  measured  data  vs  sea  model 
%  Generate  a  random  number 
U  =  rand(l); 

%  Retrieve  alpha,  beta,  and  gamma  of  Weibull  distribution  for 

%  the  requested  angle 

row  =  find(uldo(:,6)==(maxangle-l)); 

alpha  =  uldofrow,  12); 

beta  =  uldo(row,13); 

gama  =  uldo(row,14); 

%  Calculate  Sigma 
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sigma  =  beta*((-log(U)).Aalpha)+gama; 

plateau(2)  =  20*logl0(sigma)  -  seamodel(maxangle-l, sea, freq, pole); 

%  Calculate  third  point  of  measured  data  vs  sea  model 

%  Generate  a  random  number 
U  =  rand(l); 

%  Retrieve  alpha,  beta,  and  gamma  of  Weibull  distribution  for 
%  the  requested  angle 
row  =  find(uldo(:,6)==(maxangle-l)); 
alpha  =  uldo(row,12); 
beta  =  uldo(row,13); 
gama  =  uldo(row,14); 

%  Calculate  Sigma 

sigma  =  beta*((-log(U)).Aalpha)+gama; 

plateau(3)  =  20*logl0(sigma)  -  seamodel(maxangle-2, sea, freq, pole); 

%  Average  three  points  and  combine  sigmadb 
sigmadb  =  sigmadb  +  mean(plateau); 

end 

return 


A.3  SEAMODEL  Function 

function  sigmadb  =  seamodel(angle,wndsp, freq, pole) 

%  GIT  water  model  as  written  in  ALARM  91  NPS 
%  Written  23  Mar  99 

%  Capt  Bill  O'Connor  Thesis  work  for  AFRL/MBV  Lab  (Capt  Kelce  Wilson) 
%  Seamodel  will  return  a  sea  clutter  coefficient  based  on  model 
%  developed  by  Horst  at  el.  (Georgia  Tech  ,1978)  for  a  given  incidence 
%  angle  (degrees),  wind  speed  (m/s),  frequency  (Hz),  and  polarization 
%  (vertical  or  horizontal) 

%  S  =  SEAMODEL(angle,wind  speed,  frequency,  polarization) 

%  Developed  in  conjunction  with  WMOC  model  (WMOC.M) 

%  Sea  Model  works  on  grazing  angle  nomenclature  (horizon  =  0) 
angle  =  90-angle; 
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%  Convert  to  radians 
angle  =  angle  .*  (pi./180); 

%  Set  direction  of  wind  to  across  the  beam.  Curves  work  best 
%  in  this  configuration 
wnddir  =  (pi./2); 

%  Calculate  light  wavelength  &  average  water  wave  height  in  m 
lambda  =  299792500/freq; 
wave_height  =  0.348. *(wndsp./l  1).A1.8; 

if  lamda  >=  0.03  %  freq  <=  10  GHz 

%  Multipath  Interference 

rough  =  (14.4  .*  lamda  +  5.5)  .*  (angle  .*  wave_height) ./lambda; 

A_i  =  (rough.  A4)./(l  +  (rough  A4)); 

%  Interference  due  to  Wind  Speed  (or  wave  height) 

Q=  l.l.*(lambda+0.2).A(-0.4); 

A_w  =  ((1.94  .*  wndsp)./(l  +  wndsp./15))  AQ; 

%  Interference  due  to  wind  direction  ( or  wave  direction) 

A_u  =  exp(0.2.*cos(wnddir).*(l  -  2.8  .*  angle).*((lambda+0.02).A(-0.4))); 

%  Calculate  sigma 

sigmah=10.*logl0(3.9e-6.*lambda.*(angle.A0.4).*A_i.*A_u.*A_w); 

if  pole  ==  1  %  Horizontal  Polarization 
sigmadb  =  sigmah; 
else  %  Vertical  Polarization 
if  lamda  <=  0.01  %  freq  >  3  GHz 

sigmadb  =  sigmah  -  1.05.*log(wave_height+0.02)+1.09.*log(lambda)... 
+1 .27.  *log(angle+0.00 1  )+9.7 ; 
else 

sigmadb  =  sigmah  -  1.73.*log(wave_height+0.02)+3.76.*log(lambda)... 
+2.46.*log(angle+0.001)+22.2; 
end 

end 

else  %Freq  >  10GHz 

Q  =  1.93.*(lambda+  0.2).A(-0.4); 

A_w  =  ((1.94  .*  wndsp)./(l  +  wndsp./15)).AQ; 
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rough  =  (14.4  .*  lambda  +5.5)  .*  (angle  .*  wave_height)  ./lambda; 

A_i  =  (rough.A4)./(l  +  (rough.  A4)); 

A_u  =  exp(0.25.*cos(wnddir).*(l  -  2.8  .*  angle). *((lambda).A(-0.33))); 

sigmah  =  10.*logl0(5.8e-6.*(angle  A0.47).*A_i.*A_u.*A_w); 

if  pole  ==  1  %  Horizontal  Polarization 
sigmadb  =  sigmah; 
else  %  Vertical  Polarization 
sigmadb  =  sigmah  -  1.38.*log(wave_height)+3.42.*log(lambda)... 
+1.31.  *log(angle)+ 18.55; 
end 


end 

return 
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